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two  wave  fans,  each  consisting  of  rarefactions,  shocks  and  contact  discontinuities,  separated 
by  constant  states.  At  shocks  the  self-similar  viscous  solution  has  the  internal  structure  of 
traveling  waves,  and  an  admissibility  criterion  identified  by  Wendroff  [W’]  is  fullfillcd. 
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ELASTIC  AS  LIMIT  OF  VISCOELASTIC  RESPONSE, 
IN  A  CONTEXT  OF  SELF-SIMILAR  VISCOUS  LIMITS 


1.  INTRODUCTION 


The  equations  describing  one- dimensional,  isothermal  motions  of  elastic  materials  in 
Lagrangean  coordinates  are 

diu  —  dj:V  =  0 

(E) 


dtu  —  dxv  =  0 
div  —  dxcr(u)  =  0  . 


Under  the  assumption 


(Tu(ti)  >  0  (If) 

they  form  a  pair  of  strictly  hyperbolic  conservation  laws  with  characteristic  speeds  Ai(ti)  = 
—  V^cr.da),  A.;(/0  =  +  )•  As  a  model  for  nonlinear  elastic  response  (E)  provides  an 

approi)riat('  description  for  both  longitudinal  and  slu'aring  motions.  In  the  former  case  ii 
stands  for  the  longitudinal  strain  (or  specific  volume)  and  takes  strictly  positive  values, 
a{ti)  is  th(?  (ext('nsiv('  or  coini)r(.'ssive)  normal  stress  and  v  is  the  velocity.  In  the  latter 
ca.se  u  describes  the  shear  strain,  taking  now  values  on  the  whole  real  line,  (t(u)  stands  for 
the  shear  stress  and  v  for  the  velocity  in  the  shearing  direction. 

It  is  Wf.'ll  known  that,  due  to  wave  breaking  phenomena  pertaining  to  nonlinear  elastic 
response,  (E)  does  not  in  general  admit  globally  defined  smooth  solutions.  One  approach 
for  obtaining  admissible  weak  solutions  is  to  view  elastic  as  the  limiting  case  of  viscoelastic 
response  and  attempt  to  construct  solutions  of  (E)  as  limits  of  solutions  of  associated 
viscoelastic  models,  such  as 


dt «  —  5x  y  =  0 


(VE) 


dtv  -  dj;cr{u)  =  £  di{k{u)dtv)  , 

as  the  viscosity  £  -»■  0;  above  k{u)  =  1/u,  u  €  (0,oo),  for  longitudinal  motions  while 
k{u)  =  1,  u  €  (-00,00),  for  shearing  motions.  Despite  recent  advances  regarding  viscous 
limits  via  the  method  of  compensated  compactness  {c.g.  DiPerna  [Dp]),  the  understanding 


of  the  limiting  process  remains  incomplete  even  for  the  case  of  Riemann  data 


u_  ,  X  <  0 
u+  ,  X  >  0 


X  <  0 
X  >  0 


(1.1) 


r 


□  □ 


In  a  notable  exception,  for  data  corresponding  to  a  single  shock  in  (E),  HofF  and  Liu  [HLj 
obtain  detailed  information  on  the  qualitative  properties  of  the  viscous  solution  and  justify 
the  inviscid  limit. 

Due  to  the  invariance  under  dilations  of  coordinates  {x,t)  — ♦  (ax,  at),  a  >  Q,  the 
Riemann  problem  for  (E)  admits  solutions  of  the  form  (u(y),  functions  of  the 

single  variable  if  =  f.  These  are  constructed  as  weak  solutions  on  (—00,  cc)  of  the  problem 


-f 


=  0 


(1.2) 


li(±Oo)  =  ,  t'(±cc)  =  v± 


(1.3) 


subject  to  appropriate  admissibility  criteria  at  shocks  (Tupciev  [Tui],  WendrofF  [W],  Liu 
[L]  for  general  strictly  hyperbolic  systems). 

Our  objective  is  to  justify  (1.2  —  1.3)  .as  the  c  — ♦  0  limit  of  solutions  to  the  boundary 
value  problem 


(1.4) 


with  data  (1.3),  and  to  investigate  the  admissibility  restrictions  imposed  by  the  limiting 
process.  The  problem  (1.4  —  1.3)  is  denoted  by  {V.).  It  is  studied  with  the  intent  to  reveal 
aspects  of  the  more  complex  relationship  between  (VE)  and  {E). 

The  study  of  self-simil.ar  viscous  limits  was  proposed  in  the  articles  by  Kalasnikov 
[Ka]  and  Dafermos  [Di].  It  is  motivated  by  introducing  an  artificial  regularization 


di  u  —  d^v  =  0 

div  -  5xcr(u)  =  di{k{u)dtv) 


(1.5) 


that  preserves  the  invariance  under  dilations  of  coordinates.  Problem  (1.5  —  1.1)  admits 
as  solutions  functions  of  the  single  variable  f  =  x/t,  constructed  by  solving  (7^«).  The 
approach  leads  to  the  study  of  nonautonomous  boundary  value  problems  such  as  {Vc), 
and  invoh'es  variation  estimates  for  effecting  the  ;  — ►  0  limit  to  solve  the  Riemann  prob¬ 
lem.  The  procedure  has  been  carried  out  for  various  systems  of  two  conservation  laws, 
both  hyperbolic  [Di,  D2,  DDp.  KKr,  STz]  as  well  as  of  mixed  type  [Si,  Fa].  The  two 
regularizations  have  been  directly  compared  for  Burgers’  equation  [S^]. 
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We  adopt  the  following  hypotheses  on  constitutive  functions.  For  shearing  motions 
<r(u)  is  a  C^-function  defined  on  the  real  line  that  satisfies  {H)  and  has  the  behavior 

cr(u) — oo  as  u-»-oo,  cr(u)  — >  oo  zs  u -*  oo  .  {H)s 

For  longitudinzil  motions  <t(u)  is  C^-function  on  (0,  oo)  satisfying  (H)  and 

cr{u)  —*  — oo  as  u  — >  0+  ,  —^cc  as  u  —*  oo  .  {H)i 

The  latter  guarantees  that  infinite  compression  is  associated  with  infinite  compressive 
stress.  It  also  requires  that  infinite  extension  is  associated  with  infinite  extensive  stress, 
what  excludes  the  case  of  gases.  The  precise  form  of  A:(a)  has  no  influence  on  our  analysis 
as  long  as  k{u)  is  a  strictly  positive  smooth  function;  in  the  models  under  study  k{u)  =  1 
for  shearing  motions  while  k{u)  =  1/ ii  for  longitudinal  motions.  The  Riemann  data  u±.  v± 
are  taken  arbitrary  for  shearing  motions,  but  satisfy  ,  u_  >  0  for  longitudinal  motions. 

The  goals  set  in  this  article  are  (i)  to  construct  solutions  of  the  boundary  value  problem 
(Pe),  4  >  0  fi.xed,  (ii)  to  ptiss  to  the  limit  s  0  and  solve  the  Riemann  problem,  (hi)  to 
investigate  the  structure  of  the  emerging  solution.  Existence  follows  from  an  application 
of  the  Leray-Schauder  degree  theory  to  a  construction  scheme  that  captures  the  interplay 
between  the  ’’parabolic”  and  ’hyperbolic”  terms  in  (1.4).  Based  on  a-priori  estimates 
developed  in  Section  2,  the  construction  is  carried  out  in  Section  3.  Following  that  in 
Section  4,  the  passage  to  the  s  — »  0  limit  is  performed,  in  a  framework  of  uniform  in  £ 
bounds  for  the  total  variation  and  making  use  of  Helly’s  theorem.  The  outcome  of  the 
first  two  steps  can  be  summarized  as  follows:  Under  hypotheses  {H)s  or  given  any 

data  (u±,u±),  £  >  0  the  boundary  value  problem  {Vc)  admits  a  solution.  Given  a  family 
{(we,ye)}e>o  of  such  solutions  there  exists  a  subsequence  {(«e„»Vc„)},  with  £„  — +  0,  such 
that  Ue„  —*■  u,  Ve„  V  pointwise  on  the  reals.  The  function  (u,u)  is  of  bounded  variation 
and  a  weak  solution  of  the  Riemann  problem  (1.2  —  1.3). 

The  structure  of  the  limit  function  {ii,v)  is  undertaken  in  Section  5.  The  equations  of 
elcisticity  (E)  are  covered  as  a  special  case  in  Dafermos  [D2],  on  systems  of  two  conservation 
laws  with  a  full  viscosity  matrix.  Although  several  ideas  from  [D2]  are  used,  this  work  is 
based  on  a  special  property  associated  with  the  full  viscosity  matrix  which  is  not  applicable 
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here.  Instead,  we  follow  the  strategy  used  in  Tzavarsis  [Tz]  for  the  Broadwell  model. 
One  starts  with  certain  representation  formulcis  e.xpressing  the  derivatives  as 

averaging  processes  and  shows  that  there  are  finite  signed  Borel  measures  A,  i/  so  that 

(u' ^ )  — ‘  ( A,  weak—  ★  in  measures . 

The  key  ingredient  is  then  Proposition  5.1,  which  shows  that  on  the  support  of  A,  i/  a 
function  g  related  to  the  antiderivatives  of  the  eigenvalues  is  minimized.  That,  in  turn, 
yields  a  characterization  of  the  behavior  of  the  wave  speeds  Aj(u(if))  at  places  where  the 
solution  is  nonconstant.  As  it  turns  out.  (ii,v)  is  composed  of  two  wave  fans  separated  by 
a  constant  state.  Each  wave  fan  is  associated  with  one  of  the  characteristic  fields  of  {E) 
and  consists  of  either  a  single  rarefaction,  or  a  single  shock,  or  an  alternating  sequence 
of  rarefactions  and  shocks  such  that  each  shock  adjacent  to  a  rarefaction  on  one  side  is  a 
contact  on  that  side.  .At  shocks  the  Rankiue-Hugoniot  and  a  weak  form  of  the  La.\  shock 
conditions  are  satisfied. 

In  Section  6,  we  fix  a  point  ^  of  discontinuity,  with  left  and  right  limits  (u(if — ),  — )), 
(u((f4-),  u((f+)),  and  discuss  the  relation  between  self-similar  limits  and  traveling  wave 
solutions  for  the  viscous  equations.  We  refer  to  (Tu2,  Do,  Fa,  Tz]  for  investigations  of 
this  issue  in  other  contexts.  It  turns  out  that  have  the  internal  structure  of  a 

traveling  wave  solution  of  {VE),  and  that  discontinuities  emerging  via  self-similar  viscous 
limits  satisfy  for  u(^—)  <  u(^+)  and  ^  >  0 

[cr(t/)-a(u(^-))]-^Mt^-«(M]  >  0,  Ue(u(^-),u(^+)).  (1-6) 

{<  0  for  ^  <  0).  Moreover,  if  (1.6)  holds  as  a  strict  inequality  there  is  a  single  shock 
profile  connecting  u(^-)  to  u(^-l-).  On  the  other  hand,  if  (1.6)  holds  as  an  equality  at  a 
finite  number  of  points  Uj,  j  =  1,  ...,m,  there  will  be  a  chain  of  traveling  waves  connecting 
consecutively  u{(—)  to  ui,  each  of  the  points  uj  to  the  next,  and  Um  to  u(^-f-).  This 
relation  between  self-similar  limits  and  traveling  waves  was  conjectured  by  Tupciev  [Tu2]. 
Condition  (1.6)  can  be  motivated  by  an  analysis  of  traveling  waves  and  was  identified  by 
W’endroff  [W]  zis  a  criterion  to  single  out  admissible  shocks  for  the  equations  of  elasticity. 
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2.  CONNECTING  TRAJECTORIES  I  -  A  PRIORI  ESTIMATES 
Consider  the  nonlinear  boundary- value  problem  (P,) 

-^u'  -  u'  =  0 

(2.1) 

-^v'  -  a{uy  =  £{k{u)vy 

u{±oo)  =  u±  ,  v(±oo)  =  i;±  ,  (2.2) 


on  — oo  <  ^  <  oo  for  £  >  0  fixed.  Our  objective  is  to  construct  solutions  of 

(Vs),  that  is  trajectories  connecting  the  end  states  (u_,u_)  and  (u+,u+).  Both  czises  of 
longitudinal  and  shearing  motions  are  studied  simultaneously.  For  longitudinal  motions 
the  strain  u  is  expected  to  turn  out  positive,  o(ii)  is  a  C*-function  on  (O.oc)  that  satisfies 
(H)  and  {H)l,  while  k{n)  —  For  shearing  motions  u  takes  values  on  the  reals,  <t(u)  is 
a  C‘-function  on  (— oc.oo)  that  satisfies  (ff),  and  k(u)  =  1.  The  data  ti±,  v±  are 

arbitrary  for  shearing  motions,  but  satisfy  u_  >  0  for  longitudinal  motions. 

2.1  Definition  and  Regularity  Properties  for  Solutions. 

The  system  of  nonautonomous  differential  equations  (2.1)  is  singular  at  ^  =  0.  Before 
outlining  a  construction  scheme  it  is  expedient  to  clarify  the  meaning  and  to  tmalyze  the 
regularity  of  solutions.  Consider  the  weak  form  of  (2.1) 


J (Cu  +  +  j  u<pdC,  =  0 

J {(^v  -b  <7{ii))ip' d(^  ^  ^  j  j  0  ’ 

where  €  Cg(IR),  continuously  differentiable  functions  with  compact  support. 


(2.3) 

(2.4) 


Definition.  The  pair  («(O.Ks))  ^  €  ■^Lc(^)’  ^  ^  ^ 

k(u)v'  6  Lfg^(JR)  is  a  solution  of  (Vg)  if  (w,u)  satisfy  (2.3  —  2.4),  for  emy  €  C^(IR), 

emd  lim  u(^)  =  u±,  lim  v(^)  =  v±. 

±oo  ±oo 


Lemma  2.1.  Lei  (u,v)  be  a  solution  of  {Vg).  Then 

(i)  V  are  continuous  on  IR,  v'  continuous  on  El  —  {0} 

(ii)  lim  ^o^(u(0)  +  £^*(“(0)^'(0)  exist  and  are  equal 
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(iii)  The  functions  (u,u)  satisfy 

[<“(0  +  "«)]*-/  '‘(Odi  =  0  (2.5) 

[«"(«)  +  <’(>‘(0)  +  H-{u(?))v'({)]  ‘  -  /  i,{C)rfC  =  0  (2.6) 

for  any  a,6  6  IR.  In  (2.5),  (2.6)  the  appropriate  right  or  left  limits  are  understood 
whenever  a  or  b  are  equal  to  0. 


Proof.  It  follows  from  (2.3)  that  v,  6  and  thus  u,  are  continuous.  More¬ 
over.  u  €  I^g^(ni),  are  easily  seen  to  imply  ^u{^)  -+  0  as  ^  +  0. 

Fi.x  now  a,  6  €  IR.  with  a,  6  ==  0,  a  <  b  and  consider  (2.4)  with 


n 


—  oc<^<a  —  — 
n 

1 

a - ^  <  o 

n 


1 


a  <  ^  <  b 


1 


-n«-(6-h-)l 

0 


1 


b<^<b+- 

n 

6-1 —  <  ^  <  +00  . 
n 


As  i/’n  ^  Cc  (IR)  it  cannot  be  directly  used  as  a  test  function.  However,  using  appro.ximation 
by  Cl  functions,  (2.4)  can  be  easily  extended  to  hold  for  and  yields 

n  /  (Cu  -f  <7(u)  4-  £k{ii)v')d(^  (Ci^  +  <^(^0  +  ek{u)v')d^ 

Ja-^  Jb 

+  j  Vi'n<i(i  =  0  . 


Letting  n  — >  oo  shows  that  (2.6)  holds  for  a.e.  a,  6  6  IR.  Since  v,  are  continuous,  v' 
is  also  continuous  on  HI  —  {0}  and  (2.6)  holds  for  any  a  <  b  with  a  6^0.  Letting  now 
a  — ►  0—  or  6  — ►  0-f  establishes  (ii).  A  similar  argument  shows  (2.5).  | 


In  summary,  it  follows  from  (2.5)  and  (2.6)  that  «  €  C*(IR  —  {0}),  v  €  C^(IR  —  {0}) 
satisfy  the  differential  equations 


—^u'  —  u*  =  0 
—^v'  —  <7(11)'  =5  t  (^‘(u)  v’y 


(2.7) 
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for  ^  7^  0,  and  the  jump  relations 


t;(0-)  =  u(0+) 

(7(u(0-))  +  £i'(u(0-))t;'(0-)  =  <7(u(0+))  +  £k{u{0+))v' (0+) 


(2.8) 


Next  we  analyze  the  behavior  of  {u.v)  in  the  neighborhood  of  ^  =  0  and  ^  =  ±oo. 
From  (2.7)  we  obtain 


e[k[ii)v  )  + 

which  in  turn  yields 

dr  ,  f  L 

\T 

and  upon  iiitcgratiug 

(  k{u{L\4.))r'[c\  +  ) 


/V  ,  -^^u(u)  , 


V  =  0,  ^  ^0, 


1  C~  -<7>i(»(C)) 

(^k(uiQ) 


^c} 


=  0, 


=  { 


kiuiO) 


f  1  r  C--^u(n(0)  , 

d'(u(c))  ^ 


>  0 


'‘'(()  =  -U'(a  .  (f^o, 


(2.9) 


(2.10) 


(2.11) 


where  a_,  a+  are  any  cous^mits  with  a_  <  0  and  0+  >  0.  The  above  formulas  suggest 
that  the  behavior  of  e',  u'  near  ^  =  0  and  ^  =  ±oo  is  determined  by  the  behavior  of  u(^) 
as  ^  +  0  and  as  ^  ±00,  respectively. 

To  explore  this  point  further,  assume  that  u  is  bounded 


sup  |u(<f)|  <  M  (M)s 

— oo<^<oo 

in  the  case  of  shearing  motions,  or  that  u  is  bounded  from  above  and  bounded  away  from 
zero  from  below 

0  <  6  <  u(^)  <  M  ,  ^  6  (-00,00)  (M)l 

in  the  case  of  longitudinal  motions.  Then  hypothesis  (H)  and  the  form  of  k{u)  imply 

0  <  ^-0  <  kiuiO)  <  A'o  (2.12) 
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0  <  Co  <  <Tu(u(0)  <  Aq  , 


(2.13) 


where  the  constants  fco,  A'o,  do  and  --io  depend  on  <5  and  M. 

Lemma  2.2.  Let  {u,v)  be  a  solution  of  {2.7)  that  satisfies  {M)s  for  shearing  motions  or 
{M)l  for  longitudinal  motions,  so  that  (2.12)  and  (2.13)  hold.  Let  X±  —  ±.y/aQ,  A.±  = 
±\/!5o-  Then  for  0  <  ^  <  oc+  < 


X*  -o* 

1  /  ^  \-^T^ 


/{^)\<k{uia+))\v'{a^)\-(±) 

Kq  \Q(+/ 


for  A+  <  a+  <  ^ 


for  A_<a_<,^<() 


(2.14) 


ATo  X  O;  _  / 


and  for  (^  <  a_  < 


(2,10) 


(2.1-) 


Proof.  We  show  (2.14)  and  (2.15);  the  proof  of  (2.16)  and  (2.17)  is  similar. 
Fi.x  0  <  ^  <  Q+  <  A^..  Then  (2.12).  (2.13)  imply  for  0  <  C  <  a+ 


C~  -  (»((,•))  ^  C'  -  Qq 

c  -  C 


and  thus 


sxp  {  -  / 

^  C  J  e 

''  s 


1  r*  c-  -  <T.(.<(C))  ^,1 _ f  Ai  -  q5.  r*  di 


=(^) 


^  \  th’o 


Hence,  (2.14)  follows  from  (2.10). 

Let  now  0  <  .\+  <  q+  <  ^  be  fi.xed  and  ob.scrve  that  for  ^  >  Q+ 

C--M«(0)  ..  Alu 

C  "  ^  C  ~  ^ 


s 


Hence 


exp 


1  /■'  (.'-g.KO)  ..I 
«“«)) 


and  (2.15)  follows  from  (2.10).  | 

Lemma  2.2  implies  that  solutions  ((/.  {')  of  (  P,)  with  u  bounded  have  the  following 
properties: 

(i'l  |e'(^)|  —  0.  |ii'(0)|  -*  0  as  ^  ±x 

(iii  )  =  0(  |.f  1“  ),  )  =  0(  )  as  »  On.  for  some  a  >  0. 

Therefon',  e'iO-f-)  =  (’'(0  — )  =  0  and  th«'  right  and  left  limits  f/(0-i-).  u{0  —  )  exist  and  are 
finite.  By  (2.S)  and  [H], 


a(n(()-))  =  rr{n(0~))  ,  u(O-)  =  ii(04-)  ,  (2. IS) 

and  thus  u  €  C(1R)  D  C‘(IR  -  {0}),  while  t;  €  C‘(ni)  H  C2(IR  -  {0}). 

2.2  The  Construction  Sclienre. 

Our  approach  for  constructing  solutions  of  (P,)  is  to  apply  the  Leray-Schauder  de¬ 
gree  theory  (eg.  Rabinowitz  [R])  to  a  deformation  of  maps,  that  preserves  the  regularity 
structure  outlined  in  the  preceding  section.  Degree  theory  has  been  a  successful  tool  for 
establishing  connecting  trajectories  in  problems  involving  self-similar  viscous  limits  (Dafer- 
mos  [Di],  Slemrod  [S],  Slemrod  and  Tzavaras  [STz]).  Due  to  the  absence  of  diffusion  in 
(2.1)i  a  candidate  scheme  for  constructing  trajectories  should  capture  the  interplay  be¬ 
tween  the  “hyperbolic”  and  “parabolic”  effects  in  system  (2.1)  and  be  thus  different  in 
spirit  than  the  “parabolic”  based  schemes  used  in  [Di],  [S]  or  [STz]. 

Let  A’  =  C°(— oo,co)  be  the  space  of  bounded,  continuous  functions.  X  equipped 
with  the  norm  ||u||  =  sup  |n(OI  *5  a  Banach  space.  Consider  the  sets 

-oo<^<oc 

D  =  {«  €  A'  :  0  <  ^  <  uiO  <M  +  l,  -oo  <  ^  <  oo}  (2.19)/, 
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in  the  case  of  longitudinal  motions, 


Q  =  {u  ^  X  :  |u(^)l  <  M  +  1  ,  — oo  <  ^  <  00}  (2.19)s 


in  the  case  of  shearing  motions,  with  S,  M  >  0  constants  to  be  determined  later.  In  either 
case  fi  is  a  bounded  open  set  in  -Y.  For  /i  €  [0, 1],  £  >  0  fixed  and  U  €  Cl  consider 


-v'  =0 


—  00  <  if  <  oc  (2.20) 


H(±y:)  =  u±{fi)  ;=  u_  +  /t(u±  -  ii_) 


(2.21) 


c(±2c)  =  r±lji)  :=  i’_  -r  fi(r±  -  t'_)  . 

Xoro  that  for  longitudinal  morions  u  +  .  ti-  >  0  and  thus  a_(/0  >  0  for  fi  €  [0. 1]. 

Lot  ;  [0, 1]  X  Q  —  A  h('  th('  map  carrying  (fi.U)  to  the  first  component  of  the 
solution  (11 )  of  (2.20  -  2.21 ).  Then  solutions  of  (  P.. )  may  b('  visualized  as  fixed  points 
of  /"(I.-)  in  Q.  In  the  following  sc'ction  w(?  show  that  J-  is  well  dttfined  and  apply  the 
Leray-Sch.auder  degree  thetny  to  tint  homotopy  /  —  Towards  that  it  is  necessary 

to  estimate  the  solutions  of  n  —  u)  =  0,  /t  G  [0, 1|;  equivalently,  to  establish  a-priori 
estimattts  for  solutions  of  (2.1),  (2.21).  These  are  pursued  next. 


2.3.  A-priori  Estimates. 

For  the  remainder  of  the  section,  let  (a(f).t’(s))  fie  a  solution  of  (2.1)  on  (  — oc,oo), 
with  smoothness  ft  €  C(IR)  H  C^(IR  —  {0}),  v  €  C^(IR)  H  C^(IR  —  {0}).  satisfying  the 
boundary  conditions  (2.21)  for  some  /t  €  [0, 1].  For  longitudinal  motions  it  is  also  assumed 
that  u(f)  >  0  on  (— oc,oo).  In  the  sequel  C,  will  stand  for  a  generic  constant  depending 
on  Ui,  v±  and  £  but  independent  of  n,  while  C  will  denote  such  constants  whenever  they 
are  independent  of  £. 

From  (2.10)  we  see  that; 

(i)  If  y'(fo)  =  0  for  some  fo  €  (0,oo)  then  t’'(f)  =  0  on  (0,oo). 

(ii)  If  u'(fo)  =  0  for  some  fo  €  (— oo,0)  then  o'(f)  =  0  on  (— oo,0). 

Thus  on  each  of  the  segmf'iits  (— oc,0)  and  (0.  cc)  either  v  is  strictly  monotone  or  it  is  a 
constant  state.  By  (2.11)  the  same  property  characterizes  u.  Consequently  we  can  classify 
the  possible  shapes  of  («,v)  into  five  distinct  categories: 
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C\  \  u  is  increasing  on  (—00,00),  u  is  increasing  on  (— oo,0)  and  decreasing  on  (0,oo). 

C2'  u  is  decreasing  on  (—00,00),  v  is  decreasing  on  (— oo,0)  and  increasing  on  (0,oo). 
Cy.  u  is  increasing  on  (—00, 0)  and  decreasing  on  (0, 00),  v  is  increasing  on  (—00, 00). 

Cy.  u  is  decreasing  on  (— oo,0)  and  increasing  on  (0,co),  v  is  decrezising  on  (—00,00). 
C5:  (u,  v)  has  the  behavior  depicted  in  Ci  —  C4  in  one  of  the  segments  (—00, 0)  or  (0, 00) 
and  is  constant  on  the  other. 

The  classification  of  shapes  is  the  main  ingredient  leading  to  L°°-  estimates  for  solu¬ 
tions  of  the  boundary  value  problems  (2.1),  (2.21).  Such  estimates  are  obvious  for  solutions 
of  class  C5,  so  we  concentrate  on  the  other  cases. 

Lemma  2.3.  For  solations  {u,v)  of  class  Ci  or  C2: 

miu{a_,U4.}  <  u(^)  <  max{u_,U4.}  ,  (2.22) 


|f(OI<C,  (2.23) 

for  —00  <  ^  <  oc.  The  constant  C  is  independent  of  n  and  s. 

Proof.  We  present  the  proof  for  solutions  of  class  Ci,  namely  u  is  increasing  from 
u_(p)  to  u+(/i)  and  v  is  increasing  on  (— 00, 0)  and  decreasing  on  (0,oo).  Then  (2.22) 
is  clear,  and  to  show  (2.23)  it  suffices  to  bound  i;(0)  from  above. 

Since  (u,  u)  is  a  solution  of  (2.7)  satisfying  {M)s  or  >  Lemma  2.2  implies  t;*(0)  =  0 
and  t;'(0  — ♦  0  as  ^  ♦  ±00.  Integrating  (2.1)2  over  (0,oo)  we  obtain 

-  /  Cy'(C)<^C  =  <7(w+(^0)  -  <7(»(0))  <  <7(«+)  -  <^(“(0))  . 

Jo 

and  thus 

v{l)  =  v+{^i)  -  u'(C)dC  <  v+(/i)  -  ^  (iv\Qd(; 

<  max{u_,u+}  -I-  <7(u4.)  —  <7(u(0))  . 


Finally,  from  (2.1)i, 

u(0)  =  i;(l)  +  j  (,n\C,)dC^ 

Jq 

<  max{v_,  -i-  <7(u+)  —  <t(u(0))  -f  u(l)  —  u(0) 
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is  bounded  independently  of  fi  and  e.  | 

Lemma  2.4.  For  solutions  (u,u)  of  class  C3  or  C4, 

min{u_,u+}  <  u(0  <  ma.x{u_,u+}  (2.24) 

for  —00  <  ^  <  00.  Moreover: 

(i)  Under  hypotheses  (H)  and  {H)r  for  shearing  motions 

!u(OI  <  C  (2.25)s 

fir  —00  <  ^  <  00,  where  C  is  independent  of  ^  and  £. 

(ii)  Under  hypotheses  {H)  and  [H)l  for  longitudinal  motions 

Q<6<  u{i)  <  C  (2.25)l 


for  — oc  <  ^  <  cc.  with  6  and  C  positive  constants  independent  of  and  e. 
Proof,  First  we  consider  solutions  (im*)  of  class  C3.  They  satisfy  on  (— oc,oo) 

y-(/t)  <  <  v+{fi) 

min{ti_(/t),a+(/i)}  <  u(0  <  “(0) 

and,  in  view  of  (2.21),  it  suffices  to  show  that  u(0)  is  bounded  from  above. 
Integrating  (2.1)i  over  (0,oo)  we  obtain 


-  f  C“'(C)f^C  =  v+{n)  -  t;(0)  <  t;+  -  u_  . 

Jo 


In  turn,  this  yields  for  any  ^  >  0 


u(0  =  “+(/^)-y  “'(C)dC 


(2.26) 


<«+(A^)-^"c«'(C)dC  (2.27) 

<  max{u_,u+}  4-  ^  (u+  -  u-)  . 

Let  6  be  two  points  with  0  <  ^  <  ^.  Integrating  (2.1)2  over  {^,6)  leads  to 

(7(u(0)  +  £t(u(0)i-'(0  =  <^(«(^))+--^-(«(^)y(^)  +  J  Ct^'(CMC-  (2.2S) 
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Choose  Q  €  [1,2]  such  that  v'{9)  =  t;(2)  —  t;(l)  <v+—  u_.  Then 

<t(u(^))  +  £A:(u(^)y(5)  +  Cv'(CMC 

<  <t(u($))  +  £(v+  -  v-)k(u(6))  +  2(y+  -  u-) 

-  ’ 
in  conjunction  with  (2.27)  and  (2.28),  implies  that 

£  A:(u(0)  u'(0  +  <^(^(0)  <  ^  ,  for  0  <  ^  <  1  ,  (2.29) 


with  A  a  constant  independent  of  /i  and  £.  For  solutions  of  class  C3  it  is  y'  >  0  on  (O.oo) 
and  (2.29)  gives  as  ^  ►  0 

a(u(0))  <  .4  .  (2.30) 

Under  hypotheses  (H)s  or  {H)i,  (2.30)  implies  that  u(0)  is  bounded  from  above  and 
completes  the  proof  of  (2.25). 

When  (H)s  or  {H)i  are  violated  cr{u)  increases  monotonically  to  a  limiting  finite 
value  (Too  as  li  00.  If  <7oo  >  A  the  above  argument  still  leads  to  an  e-independent  bound 
for  u(0)  and  provides  (2.25)  for  restricted  classes  of  data.  However,  if  ao©  <  A  the  estimate 
disintegrates. 

Next,  we  turn  to  solutions  of  class  C^.  For  shearing  motions  the  bound  (2.25)s  follows 
from  a  similar  argument.  So  consider  the  case  of  longitudinal  motions  for  (u,y)  of  class 
C4.  Note  that 


0  <  u  :=  min{u_,  u+}  <  u±{ft)  <  max{u_,  u+}  =:  u  ,  u_  >  , 


and  that 


<  y(0  ^  v-(m) 


«(0)  <  tx(0  <  niax{u+(/i),u_(/i)} 

on  (—00,00).  It  suffices  to  show  that  u(0)  is  bounded  from  below  away  from  zero. 
Observe  first  that  from  (2.1)i 


(2.31) 


/  C«'(C)<^C  =  v(0) -v+(/^)  <*'--«+ • 

Jo 
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Choose  n  sufficiently  large  so  that  ^(u-  —  u+)  <  ju.  Then  for  ^  >  n 

=  u+(m)  - 

>  u+(fi)  -  -  f  Cu'iQdC  >  iu  . 

n  Jc  2 


(2.32) 


Ne.xt,  choose  9  S  [n,n  +  l]  such  that  v\9)  =  t’(n  +  l)  — u(n)  >  — (u-  —  y+).  Integrating 
(2.1)2  over  (C,^)  with  0<^<n<d<n4-l  and  using  (2.32),  we  obtain 

£;-(«(0)y'(0  +  o-(it(0) 

=  <t(«(«))  +  iH4S))v'{0)  +  /  (L-'iOdC 

>  (t(u(9)}  -z(t'-  -  v+)k{ii(6))  +(n  +  l)^y(^)  -  y(o) 

>  min  (<7(7)  —  £{v-  —  v+)k{q))  —  (n  +  l)(y-  —  v+)  =:  —B 


(2.33) 


with  D  a  (possibly  negative)  constant  independent  of  n  and  £.  Since  y'  <  0  on  (0,  oo), 
(2.33)  in  conjunction  with  {H)l  imply 


u(0)  >  <r”*(— S)  >  —00 , 


(2.34) 


and  (2.25)£,  follows  from  (2.31),  (2.34).  | 

The  restrictions  on  shapes  of  solutions  together  with  the  X®®-estimates  in  Lemmas 
2.3  and  2.4  imply  that  families  {(its,ye)}->o  of  solutions  to  {Vs)  are  of  uniformly  bounded 
variation. 

Corollary  2.5.  Let  {(ue,Ue)}g>o  be  a  family  of  solutions  to  the  boundary-value  problem 
(Vs)  corresponding  to  fixed  data  (u±,u±).  If  a{u)  satisfies  (H)  and  {H)s  or  {H)l,  there 
exist  constants  C  and  6  depending  on  the  data  but  independent  of  £  such  that 

on  (-00,00)  (2.35)s 


for  shearing  motions, 


0  <  ^  <  Ue(0  <  C"  .  |Ve(^)|  <  C  ,  on  (-00,00) 


(2.35)i 


14 


for  longitudinal  motions,  and 


3.  CONNECTING  TRAJECTORIES  II  -  EXISTENCE 

In  this  section  we  study  the  map  ^ :  (0, 1]  x  -♦  X  that  carries  {fi,  U)  €  [0, 1]  x  to 
the  first  component  of  the  solution  to  (2.20),  (2.21).  is  defined  by  (2.19)j[, 

for  longitudinal  motions  or  (2.19)s  for  shearing  motions;  the  constants  6  and  M  in 
the  definition  of  Q  are  chosen  as  the  minima  and  maxima,  respectively,  of  the  bounds  in 
Lemmas  2.3  and  2.4.  Our  objectives  are:  (a)  to  show  T  is  well  defined,  and  (b)  to  establish 
existence  of  solutions  for  the  equation  u  —  .^(/u  u)  =  0,  /x  6  [0, 1]. 

Let  T  :  Cl  -*  X.  S  :  Cl  X  be  the  maps  that  carry  U  €  Cl  to  the  solution  (tt.f)  = 
{T(U),  S{U))  of  the  boundary  value  problem  (2.20)  subject  to  boundary  data 

«(— oo)  =  0,  u(+oo)  =  £x+  —  £x_  ;  t’(— oo)  =  0,  v(-hoo)  =  —  i;_  . 

Then  (u_  +  v-  +  fiS(V))  satisfy  (2.20),  (2.21)  and  thus  we  may  express  the  map 

T  in  the  form 

=  u_  +  ^lT{U)  . 


3.1.  The  Linearized  Problem 

Properties  of  T,  S  derive  from  analysis  of  the  linearized  problem 


-  u'  =  0 


-^v'  -  a(0^'  = 


—  oo  <  ^  <  oo 


subject  to  data 


tx(-oo)  =  0  ,  tt(+oo)  =  u+ —  u_ 


u(-oc)  =  0  ,  u(+oo)  =  v+  —  V- 
and  with  a(4),  k{^)  continuous  functions  satisfying  on  (—00,00) 

0  <  oo  <  a(^)  <  -4o 

0  <  to  <  HO  <  Ko  . 


(3.1) 

(3.2) 

(3.3) 
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The  solution  (u,  v)  of  (3.1),  (3.2)  can  be  calculated  explicitly.  Observe  that  u'  satisfies 
the  differential  equation 


and  thus 


v'(0  = 


f  ^  J 

f  1  -  a(s)  j  ) 

c/i  sk{s)  "^"1 

1  I 

1  -  a{s)  '1 

i-iL  sLu) 

-c+i/+(^),  e>0 

s 

^<0 


with  c+,  c_  nihltrmy  coustauts. 


(3.4) 


(3.5) 


Lemma  3.1.  /+,  /_  ham:  the  followiiuj  behavior.  There  exist  constants  a,  0,  7,  8  and 
Cf  dependmg  on  n^J,  .Ao,  A-o,  A*o.  C.  aAo  depending  on  s,  such  that 


grltl-  <  4(0  <  C,|{|?  ,  0  <  10  <  1  , 
<  4(s')  <  C.e-5«’  ,  10  >  1  . 


Proof.  By  virtue  of  (3.3),  for  0  <  ^  <  1 


sk{s) 


ds 


2A’o  Aq 


while  for  ^  >  1 


7,  st(s)  7i  «(s)  J\  »*•■(«) 


(3.6) 

(3.7) 


<  -^(t"  -  1)  +  +  M 

Sl\o  Kq 

for  some  0  <  and  M  >  0.  These  bounds  together  with  (3.4)  show  one  part  of 
(3.6),  (3.7).  The  rest  follows  by  similar  arguments.  | 
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It  follows  from  (3.6),  (3.7)  that  u',  u'  are  integrable  on  (—00,00)  and  thus  (u,i;)  can 
be  calculated  by  the  formulas 


«(0  = 


c.  r 

J  —00 

-.,-)-c+  y”/ 


“(0  = 


(3.8) 


<iC  ,  «<0 


i+iOdc  ,  e>o 


f  /-(C)<<C  ,  « <  0 . 

J  —00 


(3.9) 


The  constants  c+,  c_  are  evaluated  by  requiring  continuity  of  (u.i»)  across  ^  =  0.  That 
implies 


ruio,,.  r  . 

c+  —r-di,-c-j  - —dC,  =  -{u^-u-) 

Jo  S  J-00 

C+  f  I+{0(K  +  c-  [  /-(C)dC  =  f+  -  . 

Jo  J-00 


(3.10) 


As  a  consequence  of  (3.6),  (3.7)  the  determinant  A  of  the  linear  system  (3.10)  is  bounded 
from  below  by  a  positive  constant  and  thus  (3.10)  admits  a  unique  solution  (c+,c_). 

We  turn  now  to  estimate  the  solution  (u,v)  defined  by  (3.8),  (3.9).  Let  Ke  stand  for 
a  generic  constant  depending  on  Co,  Aq,  ko,  A'o,  £  and  u±,  v±.  Then  (3.6),  (3.7),  (3.8) 
and  (3.9)  imply 

|c+I  +  (c_I<A\  (3.11) 


Therefore 


l«(OI  +  k(OI  <  (-00,00) 

MOI  +  |v'(«l  <  ,  Kl  >  1 

|t>'(OI  <  K.  ,  |u'«)|  <  /wKI?-'  ,  0  <  k|  <  1 


(3.12) 

(3.13) 

(3.14) 


with  Pj  =  oo  if  £  <  Q  =  and  pe  <  1/(1  “7)  if  Q  <  £.  We  conclude  that  v  is 
Lipschitz,  and 

l>‘«)  -  “(S)!  <  llx'lll.-.  1C  -  (3.15) 
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u  is  Holder  continuous  with  exponent  l-^<f  if  Q  <  e  and  exponent  1  (Lipschitz)  if 
e  <  a.  Note  that  the  regularity  improves  as  e  decreases. 

3.2.  Existence  of  Solutions  for 

We  state  the  main  theorem  of  this  section  establishing  existence  of  solutions  for  (Vg): 

Theorem  3.2.  Let  e  >  0  be  fixed  and  assume  that  <t(u)  satisfies  (H)  and  (H)s  in  the 
case  of  shearing  motions,  or  (H)  and  (H)l  in  the  case  of  longitudinal  motions.  Then 
there  exists  a  solution  (u((^),  of  the  boundary  value  problem  (Vg)  on  (  — cc,oo).  For 
longitudinal  motions  the  solution  satisfies  u(^)  >  0  on  (— co,oo). 

Proof.  Let  T,  S  :  Q  —*  X  be  the  maps  carrying  U  €  to  the  solution  {T{U),  S{U)) 
of  (2.20),  (3.2).  T,  5  are  properly  defined  by  (3.S).  (3.9)  and  (3.10),  where  I±  are  given 
by  (3.4)  with  k{^)  =  k{U{^)).  a{^)  =  a{U{^)).  Define  the  map  .F  :  [0, 1]  x  D  — +  A'  by 
u)  =  ii_  4-  f.iT{U).  If  u  is  a  solution  of  u  =  u_  -f  T{u)  in  Q  and  we  set  v  =  S{u)  then 
{n,v)  is  a  solution  of  {Vg)  with  u  ^  Cl.  We  will  apply  the  Leray-Schauder  degree  theory 
(Rabinowitz  [R,  Ch  V])  to  solve 

u-^r(u)  =  u_,  ^€[0,1].  (3.16) 

First  we  show  T  :  Q  X  is  compact.  This  follows  from  (i),  (ii)  below. 

(i)  T{Cl)  is  precompact  in  X. 

Consider  a  sequence  {Un}  C  Cl  and  let  u„  =  T{Un)-  Estimates  (3.12),  (3.15)  and  (3.14) 
imply  that  the  sequence  of  functions  {wn}  are  uniformly  bounded  and  uniformly  equicon- 
tinuous  on  \—L,L]  for  any  L  >  1.  Also, 

MOI<C.e-^'  ,  ^<-1 

|u„(0-(u+-u_)|<C.e-|«’  ,  i>l, 

with  Cg  independent  of  n.  It  follows  from  the  Ascoli-Arzela  theorem  and  a  diagonal 
argument  that  there  is  a  subsequence  {u,n  }  and  a  continuous  function  u  such  that  Un*  —*  u 
uniformly  on  compact  subsets  of  the  reals.  But  then  (3.17)  implies  that  u,,*  — ♦  u  in  A*. 

(ii)  T  :  Q  —*  X  is  continuous. 
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Let  {Un}  C  be  a  sequence  such  that  -*  U  in  X  and  set  u„  =  T(Un)-  Then  u„,  U„ 
are  related  by 


=  («+-“-) +  4  J 


mo 


dC 


mo  = 


«  c 

1  -  a{U„(s)) 


(3.1S) 


ds 


kauo) ^  « A 

for  ^  >  0  and  corresponding  formulas  for  ^  <  0;  c^.,  c"  solve  (3.10).  By  part  (i)  there  is 
a  subsequence  {u/^}  and  v  6  A'  such  that  -♦  u  in  A'.  Using  (3.6),  (3.7)  we  pass  to 
the  limit  in  (3.10),  (3. IS)  along  the  subsequence  {L’n*}  and  obtain  v  =  T{U).  Since  any 
limiting  point  of  {n,i}  in  A*  is  of  the  form  v  =  T{U)  we  deduce  T{Un)  —*  T{U)  in  X. 


The  map  /tT  :  [0. 1]  x  Q  — *  A*  is  compact,  thus  the  Leray-Schauder  degree  of  /  — 
is  well  defined.  By  Lemmas  2.3.  2.4  and  the  definition  of  Q  any  solution  u  of  (3.16)  lies  in 
the  interior  of  Q.  Therefore 


d(/-/tT,Q,»_)  =  d(/,fi,«_)  =  1  ,  /i  €  [0,1]  , 
and  (3.16)  admits  at  least  one  solution  for  each  jx  €  [0, 1).  | 


4.  SOLUTION  OF  THE  RIEMANN  PROBLEM 

For  the  remainder  of  the  article  we  study  the  t  — ►  0  limit  of  solutions  to  {Vc).  Let 
{( «,.  i;5))f>o  be  a  family  of  such  solutions,  corresponding  to  fi.xed  boundary  data  (u±,  u±), 
and  having  the  property 

(ii,,Ue)  satisfy  the  uniform  bounds  (2.35),  (2.36)  for  any  e>0.  (A) 

Belly’s  selection  principle  implies  there  exists  a  subsequence  of  the  original  family,  denoted 
again  by  {(us,Ue)}  with  £  — »  0,  and  a  pair  of  functions  (u,u)  of  bounded  variation  such 
that 

.  Ve(0 -♦  y(0>  pointwise  on  (-00,00)  .  (4.1) 

Our  objectives  arc  to  show  {u,v)  is  a  solution  of  the  Ricmann  problem  and  to  investigate 
its  structure. 

Concerning  Hypothesis  (A)  we  remark;  First,  Corollary  2.5  guarantees  (A)  for  families 
corresponding  to  any  data  provided  <t(u)  satisfies  (H)  and  (/f)L  or  (H)s.  Although  (A) 
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is  not  expected  to  be  in  generzil  valid  when  (if)s  or  {H)l  are  violated,  Lemmas  2.3  and 
2.4  suggest  (-4)  still  holdt  under  relaxed  conditions  on  <j{u)  provided  the  data  (u±,  u±)  are 
appropriately  restricted.  To  account  for  such  cases,  we  work  under  the  framework  of  (-4) 
for  functions  <t(u)  that  satisfy  (if).  We  refer  to  Dafermos  and  DiPerna  [DD],  Keyfitz  and 
Kranzer  [KKr]  for  studies  of  c  — ♦  0  limits  in  the  absence  of  uniform  bounds;  see  Antman 
[A]  for  a  discussion  of  the  problem  of  infinite  compression  in  a  different  context. 

Along  families  {(ue,t;e)}  satisfying  (.4),  the  estimates  (2.12),  (2.13)  also  hold  and  the 
wave  speeds  Ai(u)  =  —yJou.{u)  ,  ^^(u)  =  +  v^<7u(u)  of  the  hyperbolic  system  (£)  are 
bounded  and  separated 

A-  <  A,(»,(s^))  <  A_  <  0  <  A+  <  X2(u,(0)  <  A+  (4.2) 

by  constants  \±  =  “  i\/^  in(lcp«mdciit  of  First,  we  show  that  (a,  c)  is  a 

solution  of  the  Riemanu  problem. 

Theorem  4.1.  Suppose  cr(u)  satisfies  (H)  and  lei  {(«»,  fj)}£>o  a  family  of  solutions 
of  ifPe)  coTTCspondvig  to  data  {ti±,v±)  and  satisfying  (.4).  There  exists  a  subsequence 
{(wen,^en)}  ^  0  and  a  hounded  function  (u,u)  of  bounded  variation  such  that 

— +  u,  Vg^  — »  V  pointwise  on  the  reals.  The  limit  {u,v)  satisfies  (1.2)  m  the  sense  of 
measures  and 

(u_,u_)  ^  <  A_ 

(u(0,  via)  =  («(0),  viO))  A_  <  ^  <  A+  (4.3) 

,  (u+,t;+)  ^  >  A+ . 

Proof.  Let  {(uejUj)}  be  a  convergent  subsequence  as  in  (4.1).  Integrating  (2.1)2  over 
(0,0  and  using  Ue(0)  =  0  yields 

ekMoy.io  =  f  ".(cw  -  (^,{0 + ■’(“.(0))  -  <t(u.{0) 

Jo 

and,  in  view  of  (A), 

t(“.(«)k;(OI  <  fdtl  +  1) .  «€  (-00,00).  (4.4) 
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Lemma  2.2,  together  with  (4.4),  implies  the  following  bounds  for  (u^,t;'):  Fix  first 
any  points  a  <  0  <  6  with  [a,  6]  C  (A_,A+).  We  see  from  (2.14),  (2.16)  and  (4.4)  that 


-r-) 


0<^<6 


a<f<0 

where  d  depends  on  the  distances  \a  —  A_[,  [6  —  A+|  and  degenerates  cis  these  distances  go 
to  zero.  In  a  similar  fashion  for  any  (— oc,a]  C  (— oo,A_)  and  [6,oo)  C  (A+,oo) 

(  c r  ^ - 


’1(01  < 


where  d  now  depends  on  |a  —  or  \b  —  A+|  respectively.  Using  (4.5),  (4.6)  and  (2.11), 
we  deduce  that  lr',|,  \ii\\  — +  0  uniformly  on  compact  subsets  of  (A_,A+),  (~oc,.\_)  and 
(.\+.cc),  and  that  the  limiting  (u.  e)  satisfies  (4.3). 


Ob.scrve  lu^xt  that  if 


then  (2.1)  implies 


1  ->  r* 

+y  1  </(u,u)  = -V(7(u) 


+  4-<z(«e,f.,)  =  ev.{k{u^yj  . 
d^  Os 


Integrating  (4.S)  over  any  interval  [a,  6] 


/  77(ui,y,)c/(;  +  t  /  A.-(u5)(f')"dC 
Ja  Ja 

=  (07(“s,u«)  -g(«s,t;s)  +  £0eA-(ue)u')|^ 


and  using  (-4)  and  (4.4),  we  obtain 


if  <iC<C(Ia|  +  li|  +  l).  (4.10) 

J  a 

Let  now  9,  0  €  C|°(ni)  be  test  functions  and  suppose  that  suppv?  C  [a,  b].  We 
integrate  (2.1)  against  4'i‘P  respectively  and  use  (4.1)  and  the  estimate 

ej  lc{u,y^<^'d(;  <  £(^J  l;(us)(i;0*‘^C)  J  ^'(ue)v’'^dc) 

<£*/2C(|G|  +  |6|  +  l)(^‘v’'"dc)'^' 
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to  conclude  upon  taking  e  — ♦  0  that  (u,  u)  satisfy 


f  u(CV’)'  +  vxp'dC  =  0 

J — oo 

/  uCCs?)'  +  =  0  . 

J — oo 


(4.11) 


As  (ii,u)  are  functions  of  bounded  variation  (1.2)  is  satisfied  in  the  sense  of  distributions 
and  in  the  sense  of  meaisures.  I 


The  function 

^u(  J ),(.•(  I ))  ,  (x,0  6  (-00,  cc)  X  (0,  CX5)  (•1-12) 

is  a  weak  solution  of  the  Riemann  proljlem  for  (£■).  Indeed,  (1.1)  follows  from  (4.3),  and 
that  (4.12)  is  a  weak  solution  of  [E)  can  be  established  by  an  argument  similar  to  the 
one  leading  to  (4.11).  In  fact,  a  solution  of  the  form  (4.12)  is  a  weak  solution  of  (£")  in 
(  — cc.oo)  X  (0,oo)  if  and  only  if  (ii(s)> ‘‘(s))  satisfies  (4.11)  (Dafermos  (D;))).  Theorem  4.1 
in  conjunction  with  Corollary  2.5  lead  to  an  existence  theorem  for  the  Riemann  problem. 

Theorem  4.2.  Assume  that  cr(u)  satisfies  {H)  and  {H)i  or  {H)s-  Given  any  (u_,t;_), 
(u+,i;+)  {with  u-l  u+  >  0  in  the  case  of  longitudinal  motions),  there  exists  a  pair  of 
functions  of  bounded  variation  {u{^),v{^))  defined  on  (— oc,oo)  {with  u(^)  >  6  >  0  in  the 
case  of  longitudinal  motions),  stich  that  ^^(f),  ^(f)^  w  ®  weak  solution  of  the  Riemann 
problem  for  (£*). 

Since  (u,  v)  are  of  bounded  variation,  their  domain  can  be  decomposed  into  two  disjoint 
subsets,  (—00,00)  =  C  U  S,  such  that  on  C  the  function  (u,u)  is  continuous  while  on  S  it 
undergoes  jump  discontinuities.  S  C  [A_,  A_]U[A.j.,  A+]  is  at  most  countable  and  the  right 
and  left  limits  of  {u,v)  exist  at  any  ^  €  5.  Moreover,  u  and  v  inherit  the  monotonicity 
properties  of  Ue  and  Vg  listed  in  Gi  —  O5.  We  show  next  that  at  points  of  S  the  Rankine- 
Hugoniot  conditions  are  satisfied. 


Lemma  4.3.  At  any  ^ 

-e  [«(^+)  -  «(^-)]  -  T  <^-)]  =  0 
K^+)  -  t'(^-)I  -  K«(C+))  -  <7(«(^-)]  =  0  ■ 
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Proof.  This  is  a  consequence  of  the  fact  that  {u,v)  of  bounded  variation  satisfies  (4.11). 
We  outline  a  proof,  which  is  in  the  spirit  of  self-similar  viscous  limits. 

The  shapes  of  (ue,u<)  together  with  the  Ir®“-bounds  imply 


r  \<\d<:+  r  wjdc<c. 

•/  — OO  y  — oo 


(4.14) 


The  solutions  (u«,  Vg)  of  (Ve)  satisfy  (2.5  —  2.6).  For  ^  6  S  fixed  and  any  <5  >  0  these  imply 


6ue{6)  +  Vt(9)  d9  —  J  J  d(^d9  =  Svc{0) 

ri+i 

j  9v,{9) -r  (7iu,(9))  d9 -h  £  J  kM9))v'^{9)  d9 


(4.15) 


f(+6  j-d 

-J  j  V:{Q  d(:;d9  =  6c7{u,(0}  . 


Take  first  £  — *  0,  using  (4.1),  (.4),  (2.12)  and  (4.14),  and  then  divide  by  6  and  take  6  — ►  0- 


to  obtain 


0‘(s  +  )  +  ^s  +  )  -  /  «{CWC  =  1^(0) 

Jo 

M^+) -l-(7(«(^-l-))  -  /  i;(C)dC  =  (r(u(0))  . 

Jo 


(4.16) 


In  a  similar  manner  we  establish 


^u(^-)  4- 1;(^-)  -  f  uiOdd;  =  i;(0) 

Jo 

+  <7(u(^-))  -  /  v(CMC  =  <7(w(0))  1 
Jo 


(4.17) 


and  (4.13)  follows  from  a  comparison  of  (4.16)  and  (4.17).  | 


5.  STRUCTURE  OF  SOLUTIONS  OF  THE  RIEMANN  PROBLEM 

Let  {(ug,  Ue)}  be  a  sequence  of  solutions  of  (Pg)  as  in  Theorem  4.1  converging  pointwise 
to  the  function  (u,u)  as  £  — »  0.  We  study  next  the  structure  of  (u,u). 

First  certain  representation  formulas,  e.xpressing  the  derivatives  (u',y')  as  averaging 
processes,  are  derived.  Let 


CcU)  = 


(4  -  a„(ii,({))  _  1 

(  !:(“.(«))  «  '•■(“.(0) 


(«-A,(u,(«))(f-A,(u,(0))  (5.1) 
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= 


c,(CMC . 


[  ^  Ce(CK  ,  e  <  0 

where  a_  <  0  <  a^.  any  fixed  constants.  Integrating  (2.10)  we  obtain 


^  >0 


-v.(0)  =  k{^,{a+))v’M 

.,(0)  -  V-  =  k(u.(a.))v'.(a-)  • 

In  turn,  using  (2.10)  and  (2.11),  we  obtain  a  representation  formula  for  u' 


(tv  -t-.(O)) 


= 


_ l—e-TS.tsO 

kiti.d))  _ 

JO  A(u.(0)  ^ 


and  a  corresponding  formula  for  u'^ 


^  >0 


^<0 


<(0  = 


-Cv .i,  .  t>0 

^  J-oo 

Note  that  ge  depends  on  the  solution  (uj,Ue). 

It  is  instructive  to  use  the  correspondence  between  functions  of  bounded  variation 
and  finite  Borel  measures  on  the  reals  (  Folland  [F.  Sec  3.5]).  Let  (A,,  i/,)  be  the  measures 
associated  with  {u[,Vg)  defined  by 


measures 


<  Ae,V’  >  =  / 

\  (5.5) 

<i/e,v?>  =  y  9(^k(o<^^. 

where  0,  V’  €  Cc(ni)  continuous  functions  with  compact  support.  Then  (2.36),  (4.1)  and 
(4.14)  imply  there  exist  finite  (signed)  Borel-Stieltjes  measures  (A,t/)  such  that 

f  Ip  Ugd^  du  =:<  X,tp  >  ,  for  ^  €  Cc(IR)  , 

r  r 

/  9  Vgd^  -*  I  If  dv  =:<  >  ,  for  (j?  €  Cc(IR)  , 


(5.6) 
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that  is  Aj  — *■  A,  Vg  v  in  the  weak-*  topology  of  meeisures  (Folland  [F  ,  Sec  7.3]).  Note 
that  A  is  generated  by  the  right  continuous  function  {u{^+)  —  u_),  while  i/  is  generated 
by  It  follows  from  (2.11),  (5.3),  (5.4),  (4.5)  and  (4.6)  that  (A,^)  enjoy  the 

following  properties  ; 

(i)  supp  A  =  supp  u  C  [A+,  A_]  U  [A+,  A_]. 

(ii)  if  ^  supp  A  if  and  only  if  there  is  an  open  interval  /  3  f  such  that  u  is  constant  on  I. 

(iii)  f  ^  supp  u  if  and  only  if  there  is  an  open  interval  /  9  f  such  that  v  is  constant  on  /. 

(iv)  S  C  supp  A  =  supp  u. 

In  the  next  proposition  we  establish  an  important  property  of  the  measures  A  and  i/, 
that  incorporates  aclinissihility  restrictions  induced  by  the  self-simihir  viscosity.  To  this 
end.  it  is  convenient  to  work  with  the  quantity 


/'-<(^)  = 


Jo  t(u.(C))^ 

I.-  e-7S<iOdC 


,  f>0 


,  ^<0. 


Note  that,  since  {/ij}  are  uniformly  bounded  in  I '(El),  along  a  further  subsequence 

<^ig,ip>=  J  ;t,(f)vj(f)c/f -»< /i,V7  >,  for  V?  e  Cc(ni) ,  (5.S) 

and  that  except  for  the  trivial  cases  u+  =  t;(0)  or  t’_  =  i;(0)  it  is  supp/u  =  supp  A  =  supp  u. 
Along  the  convergent  sequence  {tig,Vg)  -*  («,y),  the  functions  Cg,  gg  converge  pointwise 
for  f  7^  0 


9:iO  — >  j(0  := 


r  c(C¥C  ■ 

Jo+ 

f  c(CR  ■ 

•/a- 


f  >0 


f  <0 


(5.10) 


We  show  that  points  in  the  support  of  f.i  are  global  minima  for  the  function  g. 
Proposition  5.1.  Let  f  G  supp/i. 

(i)  //f  €  supp/t  n  (-00,0)  then  g{0  >  g{^)  for  f  €  (-oo,0). 


(ii)  If^e  supp  ^  n  (0,  CX))  then  g{Q  >  g{0  for  C  €  (0,oo). 

Proof.  We  show  (ii);  the  proof  of  (i)  is  similar. 

Consider  the  sequence  defined  in  (5.2).  In  view  of  (5.1)  and  (2.35)  on  any  compact 
[a,  6]  C  (0,oo)  {^e}  is  uniformly  bounded  and  equicontinuous.  The  Ascoli-Arzela  theorem 
implies  a  subsequence  {ge> }  converges  uniformly.  Since  the  whole  sequence  converges 
pointwise  to  g,  we  conclude 

ge  —*  g  uniformly  on  any  [a,fe]  C  (0,oo)  .  (5-11) 

Suppose  now  that  for  some  ^  €  (0,oo)  the  conclusion  of  (ii)  is  violated.  Since  g  is 
continuous  on  (0,oo),  there  exists  some  a  >  0  such  that  the  set 

-4  ={(,'€  (0,oo)  :  <7(0 -17(0  <  -a  <  0}  (5.12) 

has  positive  Lebesgue  measure  m(>l).  We  will  show: 

Claim.  If  m(A)  >  0  there  exists  an  open  inter\'al  /  3  ^  such  that  <  >=  0  for  any 

€  Ceil). 

The  claim  implies  ^  ^  supp  /z  and  proves  (ii). 

It  remains  to  establish  the  claim.  In  view  of  (5.1),  (4.2)  and  (2.12)  there  are  positive 
constants  mi,  m2  (independent  of  s)  such  that 

C£(s)  <  for  0  <  5  <  ^ 

s  2 

Ce(s)  >  m2  (5  —  A+)  for  2A+  <  s  <  0  . 

Letting  e  — >  0  we  obtain 

,  V  rni  -  „  A+ 

c(s)  < - for  0  <  s  <  — 

S  M 

c(s)  >  m2  (5  —  A+)  for  2A+  <  s  <  0 
and,  hence,  for  0  <  C  <  ^ 

f  c{s)ds  >  —mi  — ♦  +00,  as  ^  0+ 

j 
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while  for  2A+  <  C 


[  c{s)ds  >  ^  ((C  -  A+)^  -  A+^)  +00,  as  C  +oo  . 

J2A+  ^ 

It  follows  that  points  (  near  the  boundary  of  (0,  oo)  violate  the  inequality  in  (5.12)  and 
thus  A  is  contained  in  some  compact  interval  (c,  d]  C  (0,  oo). 

Fix  now  ^  >  0  such  that 


-  <5,(f  +  ^)  implies  1^(0  -  g{^')\  <  g  • 

Since  A  C  [c,d]  and  becau.se  of  (5.11),  we  can  choose  to  >  0  such  that 

-  .5,^  +  6),  t  <  CO,  C  €  >1  imply  g'{0-9'W)<~^ 
Now  (5.7)  and  (2.12)  yield  for  I  +  S) 


0  <  fi.W)  = 


~  ih  - geW))} dc 

^  A'o 

< - e  . 

m{A) 


(5.13) 


Let  if  €  Cc(T)-  Then  (5.S)  and  (5.13)  imply 

<H.,^p>=  f  ^  0  ,  as  £ -f  0  . 

Ai-s,^+s) 

Hence,  <  /i,v?  >=  0  for  v?  €  Cc(f),  and  the  proof  of  the  claim  is  complete.  | 

The  minimization  properties  for  g  yield  information  on  the  structure  of  (u,u).  In 
particular,  a  weak  form  of  the  Lax  shock  conditions  is  induced  at  points  of  discontinuity. 
Recall  that  Ai(u)  <  0  <  A2(u)  are  the  two  wave  speeds  of  the  hyperbolic  system  (E). 

Corollary  5.2.  Let  ^,  ^'  €  HI  with  ^ 

(^)  If  ^  €  C  n  supp/<  then 


^  =  Ai(u(0)  for  ^<0, 
^  =  A2(u(0)  for  ^  >  0  . 


(5.14) 
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(b)  //^  €  <5  then  {u,v)  satisfies  at  ^  the  jump  conditions  (4.13)  and  the  inequalities 


-^i(«(^+))  <  ^  C<0, 

'^2(u(^+))  <  ^  <  >2(u(C-))  for  $  >  0  . 

(c)  If^,  €  supp/zn(-oo,0)  then  Ai(tt(^+))  =  Ai(u(^'-))  =  and  for  9  e 

e^Xxiuie))  if  See, 

\xiui9+))  =  9  =  Xiiu{d-))  if  9es. 
e  supp /tf.  ((3?^)  then  A2(u(^+))  =  A2(u(^'-))  =  and  for  9  e 

9  =  Xoiuid))  if  9eC , 

X2{u{d+))  =  9  =  X2iu{6-))  if  9eS. 


(5.15) 


(5.16) 


(5.17) 


Proof.  We  present  the  proof  for  ^  €  (O.oc).  The  function  g  in  (5.10)  is  continuous  and 
satisfies  g{^)  —*  -f  oe  as  ^  — >  0-f  or  ^  — >  +oc.  Since  c  is  of  bounded  variation 

jiin  =  litn  f  c{s)cU  =  c(^±)  ,  (5.1S) 

that  is  ^.exists  and  is  continuous  at  points  of  C,  while  only  the  right  and  left  deri\'atives 
exist  at  points  of  S.  Let  ^  €  supp/t  fl  (O.oo),  then 


<7(0  >  <7(0  for  0  <  C  <  oo 


and  thus 

c(^+)>0  ,  c(^-)<0. 

In  turn,  (5.9)  and  the  separation  properties  of  the  wave  speeds  imply  for  ^  >  0 

e-A2(«(^+))>0  ,  e-A2(u(^-))<0,  (5.19) 

which  leads  to  (5.14)  for  ^  €  C  and  to  (5.15)  for  ^  ^  S. 

It  remains  to  show'  (c).  Let  €  supp/t  H  (0,oo)  with  ^  Then  0  ^ire  both 

global  minima  for  g  with  <7(^)  =  <7(^')-  AVe  claim: 


sW  =  s(0  for  any  »€((,('). 


(5.20) 


If  (5.20)  is  violated  at  some  point,  there  exist  a,  6  with  ^  <  a  <b  <  such  that 
g{a)  =  g{b)  -  g{0  ,  g{d)  >  giO  for  a<e<b. 

At  the  points  a,  b  we  have 

A2(u(a+))  <  a  <  A2(u(a-)) 

A2(t^(6+))  <b<  A2(u(6-))  . 

On  the  other  hand  it  follows  from  Proposition  5.1  that  the  function  remains 

constant  on  the  inter\'al  (a,  6)  and  thus  A2(u(a+))  =  A2(u(6— )).  The  inequalities  then 
imply  b  <  a  which  contradicts  a  <  6;  hence  (5.20)  follows.  | 

In  summary,  the  region  where  {u,v)  is  nonconstant  consists  of  two  disjoint  closed 
intervals:  Ixi  —  [ai,6i]  C  (— oo.O)  associated  with  the  negative  characteristic  speed  Ai(u) 
and  /aj  =  [021^2]  C  (0,oo)  associated  with  the  positive  characteristic  speed  A2(u).  Each 
of  or  /a,  could  be  empty  or  consist  of  just  a  single  point.  The  function  (u,i;)  takes 
constant  values  on  the  complement  of  /a,  U  /a,  and  has  the  properties  listed  in  Corollary 
5.2  at  points  of  Ixi  or  Ix,. 

Let  /i(u,  u),  /2(u,  v)  be  the  left  eigenvectors  and  ri(u,  u),  r2(u,  v)  the  right  eigenvectors 
eissociated  with  the  eigenvalues  Ai(u),  ^ziu),  respectively,  for  the  hyperbolic  system  (E). 
They  are  normalized  to  satisfy 


li(u,v)-rj{u,v)  =  6ij  . 


The  behavior  of  the  wave  speeds  described  in  Corollary  5.2  together  with  a  proposition 
from  [Tz,  Proposition  4.6]  give. 

Proposition  5.3.  Suppose  that  /a*  =  [afc,6A.]  is  a  full  interval,  at  <  bk,  for  k  =  l  or  2. 
(i)  For  each  ^  €  [ajk,6it)  such  that  VAjt(u(^+))  •  r„(u(^+),  u(^+))  ^  0  it  is 


lim 

fc— o,A>o  h 


iWl 
1 


(5.21) 


VAA.(u(^+))T|k(ti(f+),t;(^+)) 


n(“(^+),v(^+)) 
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(ii)  For  each  ^  £  (at,6jk]  mch  that  VAi..(u(^— ))  •  rfc(i£(^— ), u(^— ))  ^  0  it  ts 


lim 

k— 0,A<0 


1 

h 


(“)«  +  A+)- 


_ 1 _ 

VAfc(ti(^-))  •  rfc(u(^-),  v{(-)) 


rk{u{^-),v{^-))  . 


(5.22) 


As  a  consequence,  (u,y)  has  right  and  left  derivatives  at  any  point  f  that  is  not 
an  accumulation  point  of  S.  If  such  a  point  ^  £  C  then  (u,u)  is  Lipshitz  there,  and  if, 
in  addition,  it  is  an  interior  point  of  /a*  then  (u.v)  is  differentiable  there.  Moreover,  a 
complete  description  of  the  structure  of  the  two  wave  fans  is  obtained.  We  distinguish  the 
following  cases: 

(i)  If  /a*  consists  of  a  single  point  then  the  solution  is  a  shock  wave  satisfying  the  weak 
form  of  the  Lax  shock  conditions. 

(ii)  If  Ix^  is  a  full  interval  of  points  in  C  the  solution  is  a  t-rarefaction  wave  (provided 
that  VAfc  •  TA-  ^  0  on  /a*  which  is  anyway  necessary  for  rarefactions). 

(iii)  In  general  /a*  consists  of  an  alternating  sequence  of  ^.--rarefactions  and  shocks  such 
that  each  shock  adjacent  to  a  rarefaction  from  one  side  is  a  contact  discontinuity  on 
that  side. 

The  emerging  picture  is  that  typical  of  strictly  hyperbolic  but  not  genuinely  nonlinear 
systems.  In  the  absence  of  genuine  nonlinearity,  the  Lax  shock  conditions  are  not  sufficient 
to  single  out  the  admissible  shocks  and  have  to  be  strengthened  (Wendroff  [W],  Liu  [L]). 
In  our  setup,  additional  restrictions  result  from  the  analysis  in  the  following  section. 


6.  SELF-SIMILAR  VISCOUS  LIMITS  AND  TRAVELING  WAVES 

In  this  section  we  discuss  the  relation  between  self-similar  viscous  limits  and  traveling 
waves  in  the  context  of  the  elasticity  equations  (E).  Let  ^  be  a  point  of  discontinuity  of 
(u,u)  and  note  that  (u(^— ),  u(^— )),  (u(^+),t;(^-b))  satisfy  the  Rankine-Hugoniot  condi¬ 
tions  (4.13).  Consider  a  sequence  of  points  with  the  property  ^  as  e  — *■  0. 

Define  functions 

c^.(0  =  “.(«.+£<).  V,(0  =  +  £  0.  (61) 
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this  introduces  a  stretching  of  the  independent  vziriable  centered  around  the  point  a 
shifted  version  of  The  uniform  estimates  (2.35),  (2.36)  imply  that  (Ue,Vt)  are  uniformly 
bounded  and  that 


TV^  U,i-)  =  TV^  u,(6  •)  ==  TV^  ue(.)  <  C  , 

TVc  K(-)  =  TV^  t;.(6  +  £  •)  =  TV^  v,{-)  <  C  . 


(6.2) 


Using  Helly’s  theorem  and  a  diagonal  argument  we  establish  the  existence  of  a  subsequence 
and  a  function  (U,  V)  such  that 


We(^«  +  £  C)  ^(C).  ^«(6  +  ^  C)  “♦  ^(C) »  pointwise  for  -  oo  <  (  <  oo  .  (6.3) 


Proposition  6.1.  Let  ^  €  S  and  suppose  that  is  a  sequence  of  points  with  .f,  — >  (f. 
Then  the  function  (C/(Q,  V’((,"))  defined  in  (6-3)  is  continuously  differentiable  and  satisfies 
on  (—00,00)  the  traveling  wave  equations 

-au-u{(-)\-[v-v(i-)\  =  o 
_4|K  _  „({_)!  -  [„^U)  -  a(«({-))|  =  HU) 

with  initial  conditions 


U{0)  =  lim  Ue(6)  >  ^(0)  =  • 


The  limits  lim  U{Q=:U±., 

C— ±00 


satisfy  the  algebraic  equations 


lim  F(C)  =:  V±  exist,  are  finite,  and  {U-,  V_),  (C’+, 
C— ±00 


14) 


-i[V-  v(e-)|  -  [HU)  -  <r(««-))l  =  0 


(6.6) 


Proof.  The  function  (ue,Ue)  satisfies  the  weak  form  (2.5),  (2.6).  First  evaluate  (2.5) 
between  the  points  and  6  to  obtain 

/•f.+eC 

(^* +£C)u«(^e+£C)  +  u*(6+£C)-^U-'W-^e(^)-y^  Us(r)dr  =  0.  (6.7) 

In  (6.7)  let  first  £  — »  0  and  then  let  consecutively  9  and  6  -+  .  Using  (6.3),  (4.1) 

and  (2.35),  we  arrive  at 

^  [uiO  -  u{^±))  +  (V(C)  -  u(^±))  =  0  .  (6.8) 
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Next,  we  evzJuate  (2.6)  at  the  points  ^  C>  ^  and  integrate  the  resulting  equation 
in  8  between  ^  and  ^  +  6,  for  some  S,  to  arrive  at 

1 

y  Ve{T)dTd6  =  -£k{Uei^e  +£(;))  V^{^e+eO 

+  £-J  k{u,i6))v'Mdd  . 


1 

+  £  CK(^«  +  C)  +  +  £  0)  -  J  y  &Ve{d)  +  <T(Ue(0))  dS 


After  some  algebraic  manipulations  and  an  integration  in  C  we  get 
(6  -i-  +  £s)) 

Jo  k{Ug{^g+£6)) 

1  /■<,'  1  /•?+^  /■s«+«a 

-si  A.))  I  I 

= + .0  -  ^.(6)  -  .(i  {[ 

Letting  e  — >  0  and  using  (6.3),  (4.1),  (2.35),  (4.14)  and  (2.12),  we  deduce 
K(5) +<r(C/»)  .  /I  . \r  ds 


f. 


k{U{s)) 
ds 


(6.9) 


From  (6.9),  by  letting  consecutively  6  — ♦  0+  and  6  — ►  0— ,  we  obtain 

I  tfC^  “  '^(“(<±)))<i“  =  -^(C)  +  ^(0) .  (6. 


10) 


It  follows  from  (6.8),  (6.10)  that  (C’(C),  V’(C))  are  continuously  differentiable  functions 
that  satisfy  the  equations  (6.4)  and  the  initial  conditions  (6.5).  Furthermore,  the  function 
U{Q  satisfies  the  differential  equation 


i k{V)  ^  -eiu-  -U-)! 

«<» 


(6.11) 
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Since  U,  V  are  of  bounded  variation  on  IR,  the  limits  lim  {U{0,  V’(C))  =:  {U±,  V±) 

C— '±00 

exist  and  are  finite.  It  is  then  straightforward  to  show  that  U+,  U-  are  equilibrixim  states 
of  (6.11) 

l.T(Cr±)  -  <r(.u(,i-))]  -eiu±-  ■<({-)  1  =  0.  (6.12) 

In  turn,  (6.8)  implies  that  (U-,  V_),  (U+,  V+)  solve  the  algebraic  equations  (6.6).  | 

The  function  (17,  V),  as  well  as  the  limiting  values  U±,  V±,  depends  on  the  choice  of 
the  sequence  For  certain  choices  of  {^4}  it  may  happen  that  {U-,  V_)  =  (i7+,  V+)  and 
the  traveling  wave  disintegrates  to  a  constant  solution.  However,  whenever  a  choice  of  shifts 
produces  a  nonconstant  solution  (U,  V)  this  will  be  a  heteroclinic  orbit  for  (6.4)  connecting 
two  distinct  equilibria  (t^_,  V'_)  and  (£■>,  V+)  solving  (6.6).  In  the  latter  case  U{C)  is  a 
solution  of  (6.11)  connecting  U-  to  what  imposes  restrictions  on  all  intermediate  states 
U.  If,  for  instance,  U-  <  U+  and  ^  >  0  then 

[,7(l/)-cr(«(?-))|-^2(i>--K(?-)]>0,  for  U  elU.,U+).  (6.13) 

We  proceed  to  investigate  the  relation  that  {(ue,t;g)}  bears  to  the  corresponding  limit 
(u,u)  at  a  point  ^  €  *5“.  Issues  of  interest  are  to  understand  what  types  of  discontinuities 
are  admissible  and  under  what  circumstances  it  is 


(V-,  V-)=  (u({-),  Kf-)).  (Cf+.  v+)  =  (t<(«+),  »(<+))  ■  (6.14) 


We  analyze  a  concrete  case;  ^  >  0,  u(^-)  <  ti(C+)  and,  in  accordance  with  the  Rankine- 
Hugoniot  conditions,  v{(—)  >  u(^+);  other  cases  can  be  treated  in  a  similar  fashion.  In 
the  case  under  study  (4.13)  and  (5.15)  take  the  form 


.2  _  -  <y(^(C-)) 

ti($+)-u(^-) 


(6.15) 


<T„(«(^+))  <  <  ‘^u(«(e-)) 


(6.16) 


They  impose  restrictions  on  the  placement  of  the  graph  of  <t(u)  relative  to  the  chord 
connecting  the  points  with  coordinates  (u(^— ),  <7’(u($— ))),  (u(^+),  a(u(^+))). 


33 


a{u) 


A-s  a  case  study  cousicU.'r  a  configuration  where  the  graph  cuts  the  chord  at  exactly  two 
intcTinediate  points  uj  and  uj  (see  Fig.  1).  Then  the  set  of  states  (f,  V)  solving  (6.C)  with 
<  U  <  consists  of  («(^-),  («i,  (^2,  ^2),  (“(^+),  v{^+))  with 

ui,  I'j  defined  by 

i  =  l,  2. 

We  will  show  that  such  a  discontinuity,  although  compatible  with  (6.15)  and  the  Lax 
conditions  (6.16),  can  not  appear  as  an  £  0  limit  of  solutions  {u,,Vg)  oi  {Vg). 

For  >  0,  u(^— )  <  u(^4-),  the  restrictions  on  the  shapes  of  solutions  dictate  that  Ug 
is  increeising  and  Vg  is  decreasing  on  (O.oo).  Fix  states  uoi,  U02  and  U03  in  the  intervals 
wi)>  («i)  ^2)  and  (u2,  n(^4-)),  respectively,  and  define  three  sequences  of  points 
such  that 

n(^ijc)  —  Uoi ,  t  =  1,  2,  3. 

According  to  Proposition  6.1,  the  corresponding  traveling  waves  Ui,  defined  by  (6.3),  satisfy 
(6.11)  with  initial  data  t’i(O)  =  uoi-  Therefore  t*i  connects  u(^— )  to  ui,  U2  connects  Uj 
to  Uf.  while  U3  connects  uz  to  u(^+).  Note  that  C^2(C)  violates  (6.13)  and  thus  such  a 
traveling  wave  cannot  appear,  except  in  the  case  ui  and  112  coincide. 
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We  conclude  that  for  a  shock  at  ^  >  0  to  arise  as  zm  £  0  limit  of  solutions  (ue,t;e) 

of  (^«),  it  must  satisfy 

[<7(::)-,rCu({-))l-{^(Cr-„«_)|  >  0,  us(u((-),u((+)).  (6.17) 

that  is  the  graph  of  cr(it)  lies  above  the  chord  joining  the  end  states.  If  (6.17)  holds  as  a 
strict  inequality  there  is  a  choice  of  producing  a  single  shock  profile  connecting  ti(^— ) 
to  u(^+).  If  (6.17)  holds  as  an  equality  at  a  finite  number  of  points  Uj,  j  =  1,  ...,m,  there 
will  be  a  chain  of  traveling  waves  connecting  u(^— )  to  ui ,  each  of  the  points  uj  to  the  ne.xt, 
and  Um  to  u((f+).  Finally,  if  (6, 17)  holds  as  an  equality  on  an  interval,  there  is  no  shock 
profile  associated  with  the  part  of  the  solution  taking  values  on  the  interval,  but  there  are 
shock  profiles  associated  with  the  complementary  part  as  before. 

The  same  argument  shows  that  when  ^  <  0  the  inequality  in  (6.17)  is  reversed,  which 
geometrically  means  that  the  graph  of  cr(w)  must  now  lie  below  the  chord.  These  conditions 
arc  sufficient  for  solving  the  Riemann  problem  by  patching  together  elementary  solutions, 
in  the  class  of  rercfaction  waves,  shocks  and  contact  discontinuities  [W,  L). 
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